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In  this  work,  a  method,  which  post-processes  the  output  of  a  control  allocation  algo¬ 
rithm,  is  developed  to  compensate  for  actuator  dynamics.  In  reaching  this  objective,  it 
is  of  great  importance  to  consider  the  real-world  applicability  of  the  approach.  In  order 
to  obtain  an  algorithm  which  can  operate  in  real-time  on  a  typical  flight  computer,  it  is 
pertinent  to  not  add  to  the  complexity  of  the  constrained  control  allocation  algorithm. 
Likewise,  it  is  necessary  to  develop  an  approach  that  is  applicable  to  both  the  saturated 
and  unsaturated  control  effector  cases.  One  method  which  can  accomplish  these  goals  is  to 
post-process  the  output  of  the  control  allocation  algorithm  to  overdrive  the  actuators  so 
that  at  the  end  of  a  sampling  interval,  the  actual  actuator  positions  are  equivalent  to  the 
desired  actuator  positions.  This  is  the  approach  taken  in  this  work  and  it  yields  a  simple, 
but  effective,  means  of  compensating  for  actuator  dynamics. 


I.  Introduction 

Numerous  control  allocation  algorithms  have  been  developed  for  aircraft  for  the  purpose  of  providing 
commands  to  suites  of  control  effectors  to  produce  desired  moments  or  accelerations.  A  number  of  approaches 
have  been  developed  that  ensure  that  the  commands  provided  to  the  effectors  are  physically  realizable.  These 
actuator  command  signals  are  feasible  in  the  sense  that  they  do  not  exceed  hardware  rate  and  position 
limits.  Buffington1  developed  a  linear  programming  based  approach  that  separately  considered  cases  where 
sufficient  control  power  was  available  to  meet  a  moment  demand  and  a  control  deficiency  case  where  the 
moment  deficiency  was  minimized.  Bodson2  developed  a  linear  programming  approach  where  the  sufficiency 
and  deficiency  branches  were  considered  simultaneously,  which  resulted  in  computational  savings.  Quadratic 
programming  approaches  have  also  been  considered.3  In  the  early  90’s,  Durham4,5  developed  a  constrained 
control  allocation  approach  called  Direct  Allocation  that  was  based  upon  geometric  concepts  of  attainable 
moment  sets.  Page  and  Steinberg6  as  well  as  Bodson2  have  presented  excellent  survey  papers  that  compare 
and  contrast  many  of  the  control  allocation  approaches  developed  over  the  last  two  decades.  Recent  work 
in  the  area  has  resulted  in  the  development  of  control  allocation  algorithms  that  can  accommodate  cases 
where  the  moments  or  accelerations  produced  by  the  control  effectors  are  nonlinear  functions  of  the  effector 
position.7,8 

A  review  of  the  constrained  control  allocation  literature  shows  that  the  coupling  effects  that  result 
from  combining  constrained  control  allocators  and  actuator  dynamics  has  largely  been  ignored,  although 
some  recent  research  has  been  performed  on  directly  including  actuator  dynamics  in  the  control  allocation 
problem.9,10  The  underlying  assumption  of  most  previous  work  is  that  actuators  respond  instantaneously 
to  commands.  This  assumption  may  at  first  seem  justified  because  in  practice,  actuator  dynamics  are 
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typically  much  faster  than  the  rigid  body  modes  that  are  to  be  controlled.  However,  interactions  between 
a  constrained  control  allocator  and  an  actuator  with  linear  dynamics  can  result  in  a  system  that  falls  well 
short  of  its  potential.  These  interactions  can  significantly  reduce  the  effective  rate  limits  of  the  system  while 
at  the  same  time  yield  control  effector  positions  which  are  different  than  the  commanded  positions.11 

The  objective  of  this  work  is  to  develop  a  method,  which  post-processes  the  output  of  a  control  allocation 
algorithm,  to  compensate  for  actuator  dynamics.  In  reaching  this  objective,  it  is  of  great  importance  to 
consider  the  real-world  applicability  of  the  approach.  In  order  to  obtain  an  algorithm  which  can  operate  in 
real-time  on  a  typical  flight  computer,  it  is  pertinent  to  not  add  to  the  complexity  of  the  constrained  control 
allocation  algorithm.  Although  the  results  in  Harkegarcl9  and  Venkataraman,  et.  al.10  overcome  some  of 
the  issues  raised  here,  it  is  not  certain  that  either  method  can  guarantee  convergence  in  one  flight  control 
system  update  for  both  the  saturated  and  unsaturated  cases.  One  method  which  can  accomplish  the  goal  is 
to  post-process  the  output  of  the  control  allocation  algorithm  to  overdrive  the  actuators  so  that  at  the  end 
of  a  sampling  interval,  the  actual  actuator  positions  are  equivalent  to  the  desired  actuator  positions.  This 
is  the  approach  taken  in  this  work  and  it  yields  a  simple,  but  effective,  means  of  compensating  for  actuator 
dynamics. 

Bolling11  has  shown  that  the  interaction  between  first-order  actuator  dynamics  and  constrained  control 
allocation  algorithms  can  be  eliminated  by  overdriving  the  actuators.  In  this  work,  the  details  of  the 
interaction  between  constrained  control  allocators  and  first-order  actuator  dynamics  are  provided.  This 
scheme  is  then  extended  to  a  second-order  actuator  dynamics  case.  Simulation  results  are  presented  for  a 
four  control  surface  vehicle  using  a  linear  programming  based  control  allocation  algorithm  that  takes  into 
account  control  effector  position  and  rate  limits.  The  results  show  that  the  desired  control  effector  positions 
can  be  obtained  by  post-processing  the  control  allocation  commands. 


II.  Effective  Rate  Limits  and  Attenuation  of  Zero-Order-Hold  Inputs  For 

First-Order  Actuator  Dynamics 


Figure  1  shows  the  system  to  be  analyzed  in  this  work.  Inputs  to  the  control  allocation  algorithm 
consist  of  a  vector  of  desired  moment  or  acceleration  commands,  ddes  £  Mn,  and  a  vector  containing  the 
current  control  surface  deflections,  8  £  Rm.  The  output  of  the  control  allocator  is  the  commanded  control 
surface  deflection  vector,  8cmd  £  Rm.  The  actuator  dynamics  respond  to  8crn(/.  to  produce  the  actual  control 
deflections,  8.  The  individual  actuators  are  assumed  to  have  uncoupled  dynamics,  hardware  rate  limits 
±<5maxi  and  position  limits  8max.  In  most  control  allocator  implementations,  rate  limits  are  taken 
into  account  by  converting  them  into  effective  position  limits  at  the  end  of  the  next  sampling  period  and 
constraining  the  effector  commands  to  respect  the  most  restrictive  of  the  rate  or  position  limits,  i.e., 

8  =  min (5 max,  8  T  8maxcAAt)  :  8_  =  max(<5min,  8  8maxcA/\t)  (1) 


where  <5  is  the  current  location  of  the  control  effectors,  8,8  are  the  most  restrictive  upper  and  lower  bounds 
on  the  effectors,  respectively,  and  At  is  the  sampling  period  of  the  digital  flight  control  system. 

In  typical  implementations,  each  element  of  the  vector  of  rate  limits  provided  to  the  control  allocation 
algorithm,  dmaXcA,  is  taken  to  be  the  true  hardware  rate  limit  of  the  corresponding  actuator.  Hence,  the 
software  rate  limit  is  equal  to  the  hardware  rate  limit  or  <5maxcA  =  8luax~  however,  as  was  shown  by  Bolling,11 
the  effective  rate  limit  of  a  scalar  system  composed  of  a  constrained  control  allocation  algorithm  and  first- 
order  actuator  dynamics  of  the  form 

j(f)  =  a 

8cmd{s)  s  +  a 


becomes 


8 


maxEFF 


=  TS 


maxcA 


(3) 


where  T  =  1  —  e~aAt,  5maxEFF  is  the  effective  rate  limit,  and  <jmaxCA  is  the  rate  limit  provided  to  the  control 
allocation  algorithm  by  the  user.  For  example,  when  a  =  20  5maxcA  =  60°^’  and  the  flight  control 
system  operates  at  50  Hz  (At  =  0.02  sec),  the  effective  rate  limit  is  19.78°  y^.  This  result  is  well  short  of 
the  ideal  60°  —  value. 
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A  few  comments  are  in  order  regarding  Figure  1  and  the  analysis  described  above.  First,  note  that 
the  instantaneous  position  limit  given  by  Equation  1  makes  use  of  a  sampled  vector  of  actuator  position 
measurements  to  compute  the  maximum  distance  that  the  actuator  can  move  in  the  next  time  instant. 
This  is  in  contrast  to  using  the  previous  value  of  actuator  command  vector,  8cmd ,  as  is  often  done  in 
simulation.  The  motivation  for  using  actuator  measurements  is  that  when  actuator  dynamics,  disturbances, 
and  uncertainties  are  taken  into  account,  the  actuator  command  vector  and  the  true  actuator  positions  will 
differ.  This  difference  can  cause  a  control  allocator  to  generate  inappropriate  actuator  commands  that  do 
not  deliver  the  desired  moments  or  accelerations.  Thus,  feeding  the  measured  actuator  position  vector  to 
the  control  allocator  has  the  advantage  of  reducing  uncertainty  in  the  actuator  position.  It  is  important; 
however,  to  be  aware  of  the  effects  of  using  the  measured  actuator  position  vector. 

In  the  same  manner  in  which  rate  limits  are  effectively  attenuated  by  actuator  dynamics,  so  are  the 
magnitudes  of  the  commanded  changes  to  control  effector  positions.  Referring  to  Figure  1,  the  desired 
situation  would  be  for  6  =  8cmd .  However,  actuator  dynamics  alter  the  command  signals  so  that,  in  general, 
8  yf  8cmd.  For  actuators  with  high  band  widths  relative  to  the  rigid  body  modes,  this  is  not  a  serious  concern. 
However,  situations  exist  where  the  actuator  dynamics  are  not  sufficiently  fast  and  need  to  be  taken  into 
account.  In  this  section,  the  effects  of  first-order  actuator  dynamics  on  the  system  shown  in  Figure  1  will  be 
discussed.  In  order  to  illustrate  the  idea  for  post-processing  control  allocator  outputs,  the  method  proposed 
by  Bolling11  that  compensates  for  both  magnitude  and  rate  limit  attenuation  will  be  described. 

The  discrete  time  representation  of  the  first-order  actuator  dynamics  equation  is  given  by 

8(tk+ i)  =  $8(tk)  +  T  8cmd(tk)  (4) 

where  <f>  =  e_aAt,  T  =  1  — e_aAt,  and  it  has  been  assumed  that  the  input  to  the  actuator  dynamics,  8cmd(tk), 
is  held  constant  over  each  sampling  period.  The  command  to  the  actuator  can  be  written  as 

8cmd(^k)  A 8crndcA{t]i )  T  S(tk)  (5) 

where  the  commanded  incremental  change  in  actuator  position  over  one  timestep  is  defined  by  A 8cmdcA  ( tk )  = 
8 cmdc  a  (^fc)  ~  8(tk)  and  where  ScmdcA(tk)  is  the  actuator  position  command  from  the  control  allocator. 
Since  the  effector  commands  are  held  constant  over  each  sampling  period,  A ScmdcA(tk)  will  appear  to  the 
actuators  to  be  a  step  command  from  the  measured  position.  Recall  that  ScmdcA(tk)  is  calculated  by  the 
control  allocator  based  upon  the  assumption  that  the  effector  will  respond  instantaneously  to  commands. 
Substituting  Equation  5  into  Equation  4  yields 

8(tk+ i)  =  <M(4)  +  T  [A8cmdcA  (tk)  +  £(tfc)]  (6) 

Since  T  <  1,  the  incremental  command  signal  from  the  control  allocation  algorithm,  A ScmdcA(tk)  is  attenu¬ 
ated  by  the  actuator  dynamics,  thus  8(tk+ 1)  y^  8cmdcA(tk)-  The  objective  is  to  find  a  gain,  M,  that  modifies 
the  output  of  the  control  allocation  algorithm  such  that  8(tk+ i)  =  ScmdcA(tk)  =  A8cmdCA(tk)  +  8(tk)-  Hence, 

8(tk+ 1)  =  +  r  [MA5cmdcA  (tk)  +  S(tk )]  (7) 

and,  solving  for  M  yields 

M  =  l  (8) 

Thus  the  actuator  command  signal  must  be  modified  such  that 

8cmd(tk)  ~  81  A8crndcA  (tk)  T  8(t-k)  =  A ScmdcA(tk)  T  S(tk)  (9) 

1  s 

where  Ts  =  1  —  e-a»*'>™At  is  the  software  scaling  factor  and  anom  is  the  nominal  bandwidth  of  the  first-order 
actuators.  Replacing  8cmd(tk)  in  Equation  4  with  8cmd(tk)  from  Equation  9  yields 

8(tk+i)  =  *&8(tk)  +  ^8cmd(tk)  =  <$>8(tk)  +  r  —  AScmdcA(tk)  +  S(tk)  (10) 

.  s 
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which,  when  a  =  anom  and  thus  T  =  Ts,  results  in  S(tk+ 1)  =  8cmdCA(tk )•  The  distinction  between  T  and 
Ts  is  made  because  it  is  possible  for  the  actuators  to  at  some  point  have  a  bandwidth  which  is  less  than 
the  nominal  value,  due  to  power  loss,  partial  failure,  etc.  In  other  words,  the  potential  exists  for  a  <  anom. 
In  this  analysis,  it  has  been  assumed  that  a„om  is  the  nominal  bandwidth  of  the  actuators  and  is  an  upper 
bound  on  bandwidth.  That  is,  if  a  /  anom,  then  a  <  anom.  In  other  words,  the  bandwidth  of  the  actuator 
(a)  cannot  be  larger  than  the  nominal  bandwidth  ( anom ).  Under  this  assumption,  T  <  Ts. 

Since  Ts  can  be  computed  from  the  known  quantities  anom  and  At,  one  can  compensate  for  command 
increment  attenuation  using  Equation  9.  For  a  bank  of  decoupled  first-order  actuators  with  nominal  band- 
widths  of  anomi,  corresponding  values  of  TSi  can  be  computed  using  Ts.  =  (1  —  e~anomiAt).  The  command 
increment  compensation  can  then  be  implemented  in  discrete  time  as  shown  in  the  block  diagram  of  Figure  2. 
In  Figure  2,  the  dashed  lines  are  for  compensation  for  second-order  actuator  dynamics,  as  will  be  discussed  in 
the  next  section,  and  therefore,  are  not  part  of  the  first-order  actuator  dynamics  compensation  scheme.  Note 
that,  for  multiple  actuators,  M  in  Figure  2  is  a  diagonal  matrix  with  the  entries  on  the  main  diagonal  being 
rsi ,  rs2 ,  •  •  •  ,  TSm ,  where  the  subscript  m  is  defined  as  the  number  of  control  effectors.  Hence,  the  magnitude 
of  the  control  allocation  command  increment  is  modified  to  counteract  the  attenuation  that  results  from 
the  interaction  between  first-order  actuator  dynamics  and  the  control  allocator.  Note  that  the  allocator 
continues  to  operate  under  the  assumption  that  the  actuator  will  respond  instantaneously  and  ensures  that 
the  commands  obey  rate  and  position  limits.  It  is  after  the  allocator  computes  a  new  set  of  commands  that 
the  command  increment  is  scaled  and  added  to  the  measured  actuator  position.  A  beneficial  consequence 
of  this  is  that  the  effective  rate  limit  will  be  equal  to  both  the  hardware  and  software  rate  limit  when  the 
magnitude  increment  is  scaled  in  this  fashion.  Thus,  when  using  this  technique,  there  is  no  need  to  adjust 
the  software  rate  limit. 

It  should  be  pointed  out  that  this  technique  effectively  modifies  the  gain  of  the  inner- loop.  In  the  typical 
situation,  where  actuator  dynamics  are  assumed  to  be  much  faster  than  the  rigid  body  modes  to  control  and 
are  hence  ignored,  a  stability  analysis  would  use  8Cmd  as  the  input  to  the  plant.  In  this  case,  it  can  be  shown 
that  for  r  <  Ts,  even  after  this  compensation  scheme  is  applied,  <5  <  8cmd .  Therefore,  something  less  than 
or  equal  to  Scmd  is  applied  to  the  plant  and  the  loop  gain  is  reduced.  The  performance  of  the  entire  control 
system  will  be  burdened  with  the  task  of  mitigating  the  effects  of  S  ^  5Cmd-  Care  must  be  taken  in  applying 
this  approach  because  it  is  possible  for  the  loop  gain  to  be  reduced. 


III.  Attenuation  of  Zero-Order-Hold  Inputs  For  Second-Order  Actuator 

Dynamics 


With  the  knowledge  gained  from  the  first-order  actuator  dynamics  case,  it  is  now  possible  to  extend  the 
previous  results  to  second-order  actuator  dynamics.  Let  the  second-order  actuator  dynamics  be  represented 
by  the  general  equation 

M 

8cmd(s) 

which,  in  state-space  form,  becomes 


2(u’rl 


(11) 


"  m ' 

0  1 

'  S(t)  ' 

+ 

0 

S^Jt)  =  A 

'  m ' 

.  <*(*) . 

_  -ul  -2 On  _ 

.  ^)  . 

k 

.  ^(*) . 

8(t)  = 

1  0 

"  m ' 

=  C 

"  m ' 

S(t) 

S(t) 

T  T58cmd(t) 


(12) 


(13) 


Using  Equations  12  and  13,  the  discrete-time  solution  to  the  second-order  actuator  dynamics  differential 
equation  becomes 


<5(ifc+i) 


$1,1  $1,2 


8(tk) 

8(tk) 


8cmd{th) 


/'*$!. 2  (/-A-  •  1 


r)dr 


(14) 
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where 


c[>  _  gA(tfc  +  l— tfc)  _ 


$1,1 

$2,1 


$1,: 

$ 


2,2 


(15) 


is  the  state  transition  matrix.  Since  A  is  constant,  the  state  transition  matrix  depends  on  the  time  difference, 
At  =  tk+i  —  tk,  and  the  explicit  dependence  on  time  has  been  eliminated.  Performing  the  operations  required 
in  Equations  14  and  15  produces 


8(tk+i)  —  Ci5(tk)  +  C'2<5(tfc)  +  C35crnd(tk) 


(16) 


where 


y-i  w n  t  ■ 

Ci  =  — -e  sin 

Ud 


ajdAt  +  arctan 


7  At 


Co  = 


OJd 


-sin 


(17) 


C3  =  — 

OJd 


ud  ■ 


7  At 


asin(u>dAt)  —  ajdCos(uidAt)\ 


U)„ 


LOd  —  Lon\/\  —  £2  and  er  =  —(u>n.  Substituting  Equation  5  into  Equation  16  produces 

d(tk+ i)  =  CiS(tk)  +  C2S(tk)  +  C3  [A5cmdCA  (tfc)  + 


(18) 


The  objective  is  to  find  a  gain  M  that  will  modify  A 5cmdCA(tk)  in  such  a  way  that  S(tk+i)  =  ScmdCA{tk)- 
Hence,  it  is  desired  to  find  M  such  that 


A 5cmdCA{tk)  +  S(tk)  —  CiS(tk)  +  C2<5(tk)  +  C3  [M AScmdCA(tk)  +  5(tfc)] 


Solving  for  M  gives 


M  = 


A5cmdCA(tk)  +  (1  —  C3  —  Ci)6(tk)  —  C2S(tk) 


(19) 


(20) 


C3ASc7ndcA  (^fc) 

Clearly,  unlike  the  simple  first-order  case,  this  compensation  is  much  more  complex.  In  fact,  this  requires 
not  only  the  position  of  the  control  effector,  but  also  the  rate  of  change  of  the  control  effector.  As  with  the 
first-order  case,  this  method  post-processes  the  control  effector  commands,  without  modifying  the  control 
allocation  algorithm.  Referring  to  Figure  2,  for  a  system  with  more  than  one  actuator,  M  would  be  a 
diagonal  matrix  with  entries  along  the  main  diagonal  being  Ml5  M2,  ■  ■  ■  ,  Mm.  Here,  Mj  would  be  computed 
using  Equation  20  and  ojrii ,  Ci  corresponding  to  the  ith  actuator.  Also  note  that  a  combination  of  first-  and 
second-order  actuators  could  be  treated  by  using  Equation  8  for  the  first-order  actuators  and  Equation  20 
for  the  second-order  actuators.  For  second-order  actuators,  the  dashed  lines  in  Figure  2  would  be  used,  due 
to  the  fact  that  the  scaling  factor,  Equation  20,  depends  on  the  actuator’s  past  position  and  rate. 

Recall,  in  the  first-order  actuator  dynamics  case,  that  due  to  a  power  loss,  partial  failure,  etc.,  it  is 
possible  for  T  <  Ts,  that  is,  for  the  bandwidth  of  the  physical  actuator  to  be  less  than  the  bandwidth  used 
to  determine  the  scaling  factors  Mj.  In  the  second-order  actuator  dynamics  case,  two  parameters  can  vary, 
namely,  the  damping  ratio  (£)  and  the  natural  frequency  (<on).  If  the  control  system  timestep  is  large  enough, 
this  could  have  an  adverse  impact  on  the  performance  of  the  system.  However,  for  a  50  Hz  update  rate,  the 
timestep  is  0.02  sec  and  the  impact  of  differences  between  the  true  C,lou  and  those  used  in  computing  Mj 
are  minimal.  Basically,  if  the  sampling  is  fast  enough,  then  the  differences  between  two  actuator  positions 
at  the  end  of  an  update  cycle  are  small  even  when  the  actuators  have  slightly  different  damping  factors  and 
natural  frequencies. 


IV.  Simulation  Results 

In  this  section,  results  from  a  simulation  of  the  system  displayed  in  Figure  2  will  be  shown.  A  rate  and 
position  constrained  linear  programming  based  control  allocator  will  be  utilized  in  this  work.  In  this  case, 
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the  control  allocation  algorithm’s  objective,  referring  to  Figure  1,  is  to  find  8cmd  such  that 

drfes  =  B5cmf/  (21) 

where  B  is  the  control  effectiveness  matrix  and  d,^  is  typically  a  set  of  moment  or  acceleration  commands 
for  the  roll,  pitch,  and  yaw  axes.  Although,  if  feasible,  the  control  allocator  will  be  able  to  find  a  8cmd  such 
that  Equation  21  holds,  the  real  test  is  to  determine  what  happens  after  the  actuator  dynamics  operate  on 
8Cmd-  Hence,  the  overall  system  goal  is  to  achieve  such  that 

ddes  =  B  8  (22) 


Equation  22  is  the  metric  upon  which  the  quality  of  results  will  be  judged.  In  this  example,  four  control 
effectors  are  present  and  the  control  effectiveness  matrix  is  fixed  at 


-0.4 

0.4 

-0.1 

0.1 

B  = 

-0.1 

-0.1 

-0.6 

-0.6 

-0.1 

0.1 

-0.1 

0.1 

where  the  elements  of  B  have  units  of  (rad / sec2) / deg .  Since  there  are  more  control  effectors  (4)  than  axes 
to  control  (3),  a  control  mixed  or  allocator  must  be  used.  In  this  work,  a  flight-tested  linear  programming 
based  control  allocation  algorithm  with  rate  and  position  limiting  will  be  used.2-7  It  is  important  to  note, 
however,  that  any  control  allocation  method  could  be  used. 

In  the  following  simulations,  only  second-order  actuator  dynamics  will  be  used.  Each  of  the  four  actuators 
have  the  following  model 


M  -  1 

8Cmd(s )  S2  +  7.071s  +  25 


(24) 


which  gives  £  =  cun  =  5,  and  k  =  1  .  Each  actuator  is  rate  and  position  limited  by  the  following  values 


8min 


8  max 


-1.5  -1.5  -1.5  -1.5  |  (deg) 
0.3  1.5  1.5  1.5  (deg) 

,  /deg' 

10  10  10  1.5  — 

sec 


(25) 


These  limits  were  selected  so  that  at  least  one  position  and  one  rate  limit  were  in  effect  at  some  time  during 
the  simulation.  This  was  done  to  show  that  the  method  developed  in  this  work  is  applicable  when  control 
effectors  are  saturated.  As  will  be  shown,  when  compensation  is  used,  actuator  1  becomes  upper  position 
limited  and  actuator  4  becomes  rate  limited  as  the  frequency  of  the  commands  increases.  The  command 
signals,  ddes >  consist  of  chirps  of  magnitude  0.025,0.075,  and  0.05  ^4  in  the  roll,  pitch,  and  yaw  channels 
and  where  the  frequency  ranged  from  0.5  —  2  Hz  as  a  linear  function  of  time  over  a  10  sec  time  interval. 

Simulation  runs  were  performed  with  and  without  compensation  for  magnitude  attenuation  due  to  second- 
order  actuator  dynamics.  Ideal  conditions  are  when  the  actuator  dynamics  can  be  represented  by  s  =  1- 
If  sufficient  control  authority  and  ideal  conditions  exist,  then  the  control  system  would  achieve  ddes  —  B<5, 
which  is  the  best  possible  performance.  Figures  3  and  4  show  d<jes  and  B<5  without  and  with  the  magnitude 
compensation  described  in  Equation  20  applied.  Clearly,  when  the  magnitude  compensation  is  not  used 
(Figure  3),  d des  4  B<5,  and  a  large  error  exists  between  these  two  quantities.  When  magnitude  compensation 
is  used  (Figure  4),  however,  ddes  =  B<5  and  near  ideal  performance  is  achieved. 

Figures  5  and  6  show  the  control  effector  commands,  8cmdCA,  and  actual  deflections,  <5,  when  magnitude 
compensation  is  not  utilized.  As  expected,  the  actual  deflections  are  considerably  different  than  the  com¬ 
manded  deflections  and  this  directly  translates  to  the  deviations  between  d^es  and  B<5  as  shown  in  Figure  3. 
On  the  other  hand,  Figures  7  and  8  display  the  control  effector  commands,  8cmdCA,  and  actual  deflections,  8 , 
when  magnitude  compensation  is  used.  In  this  case,  the  actual  deflections  are  nearly  equal  to  the  commands 
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and  as  a  result,  ddes  =  B<5.  Notice  that  control  effector  1  is  upper  position  limited  and  control  effector  4  is 
rate  limited  for  the  last  few  seconds  of  the  simulation  run.  Thus,  using  the  simple  gain  adjustment  described 
in  Equation  20  results  in  the  actual  control  deflections  being  equal  to  the  commanded  control  deflections.  It 
is  apparent  that  adjusting  the  control  effector  command  increments  can  help  to  mitigate  adverse  interactions 
between  discrete  time  implementations  of  control  allocation  algorithms  and  actuator  dynamics. 


Figure  1.  Control  allocator  and  actuator  interconnection. 


V.  Conclusions 

Interactions  between  constrained  control  allocation  algorithms  and  the  dynamics  of  actuators  can  result 
in  degraded  performance  if  not  carefully  implemented.  The  effective  rate  limit  of  a  constrained  control 
allocator  operating  in  conjunction  with  a  rate  and  position  limited  first-order  actuator  will  always  be  less 
than  the  software  rate  limit  enforced  by  the  control  allocator.  Likewise,  it  was  shown  how  the  control 
effector  commands,  from  a  control  allocation  algorithm,  are  attenuated  by  actuator  dynamics.  One  method, 
which  can  be  used  to  extract  maximum  performance  from  such  a  system,  is  to  modify  the  control  effector 
commands.  The  gains  used  to  modify  the  commands  were  computed  for  both  first-order  actuators  and  for 
second-order  actuators.  Simulation  results  show  that  significant  improvements  can  be  achieved  by  using  this 
method  to  post-process  the  commands  from  the  control  allocation  algorithm.  Benefits  of  this  method  are 
that  the  control  allocation  algorithm  need  not  be  modified,  it  requires  minimal  additional  computations,  and 
it  is  valid  for  saturated  and  unsaturated  control  effectors.  Because  the  additional  required  computations  are 
minimal,  this  method  would  be  suitable  for  implementation  in  a  real-time  environment. 
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Figure  5.  Control  Effector  1  and  2  Positions  -  Compensation  OFF. 
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Figure  6.  Control  Effector  3  and  4  Positions  -  Compensation  OFF. 
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Figure  7.  Control  Effector  1  and  2  Positions  -  Compensation  ON. 
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